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Abstract 

We consider the effect of an external bias voltage and the spatial 
variation of the surface potential, on the damping of cantilever vibra- 
tions. The electrostatic friction is due to energy losses in the sample 
created by the electromagnetic field from the oscillating charges in- 
duced on the surface of the tip by the bias voltage and spatial variation 
of the surface potential. A similar effect arises when the tip is oscillat- 
ing in the electrostatic field created by charged defects in a dielectric 
substrate. The electrostatic friction is compared with the van der 
Waals friction originating from the fluctuating electromagnetic field 
due to quantum and thermal fluctuation of the current density inside 
the bodies. We show that the electrostatic and van der Waals fric- 
tion can be greatly enhanced if on the surfaces of the sample and the 
tip there are two-dimension (2D) systems, e.g. a 2D-electron system 
or incommensurate layers of adsorbed ions exhibiting acoustic vibra- 
tions. We show that the damping of the cantilever vibrations due to 



the electrostatic friction may be of similar magnitude as the damping 
observed in recent experiments of Stipe et al [B.C. Stipe, H.J.Mamin, 
T.D.Stowe, T.W.Kenny, and D.Rugar, Phys.Rev. Lett. 87, 0982001]. 
We also show that at short separation the van der Waals friction may 
be large enough to be measured experimentally. 

1 Introduction 

A great deal of attention has been devoted to non-contact friction between 
an atomic force microscope tip and a substrate [H 121 EH IH IS] • This problem 
is related to the role of non-contact friction for ultrasensitive force detec- 
tion experiments. The ability to detect small forces is inextricably linked to 
friction via the fluctuation-dissipation theorem. According to this theorem, 
the random force that make a small particle jitter would also cause friction 
if the particle were dragged through the medium. For example, the detec- 
tion of single spins by magnetic resonance force microscopy (Bj, which has 
been proposed for three-dimensional atomic imaging T\ and quantum com- 
putation jSj, will require force fluctuations (and consequently the friction) 
to be reduced to unprecedented levels. In addition, the search for quantum 
gravitation effects at short length scale [5], and future measurements of the 
dynamical Casimir forces ^Uj, m ay eventually be limited by non-contact fric- 
tion effects. Non-contact friction is also responsible for the frictional drag 
force between two-dimensional (2D) quantum wells [TT} IT2*1 IT3]. 

In non-contact friction the bodies are separated by a potential barrier 
thick enough to prevent electrons or other particles with a finite rest mass 
from tunneling across it, but allowing interaction via the long-range electro- 
magnetic field, which is always present in the gap between bodies and can 
have different origin. The presence of an inhomogeneous tip-sample electric 
fields is difficult to avoid, even under the best experimental conditions [3]. 
For example, even if both the tip and the sample were metallic single crys- 
tals, the tip would still have corners, and more than one crystallographic 
plane exposed. The presence of atomic steps, adsorbates, and other defects 
will also contribute to the spatial variation of the surface potential. This is 
referred to as "patch effect" . The surface potential can also be easily changed 
by applying a voltage between the tip and the sample. An inhomogeneous 
electric field can also be created by charged defects embedded in a dielec- 
tric sample. The relative motion of the charged bodies will produce friction 
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which will be denoted as the electrostatic friction. 

The electromagnetic field can also be created by the fluctuating current 
density, due to thermal and quantum fluctuations inside the solids. This 
fluctuating electromagnetic field gives rise to the well-known long-range at- 
tractive van der Waals interaction between two bodies [Tl|, and is responsible 
for radiative heat transfer. If the bodies are in relative motion, the same fluc- 
tuating electromagnetic field will give rise to a friction which is frequently 
named as the van der Waals friction. 

Recently Stipe et. al. |3] observed non-contact friction between a gold sur- 
face and a gold-coated cantilever as a function of tip-sample spacing d, tem- 
perature T, and bias voltage V. The friction force F acting on the tip 
was found to be proportional to the velocity v , F = Tv. For vibration of 
the tip parallel to the surface they found T(d) = a(T)(V 2 + V^)/d n , where 
n = 1.3 ± 0.2, and V ~ 0.2 V. At 295K, for the spacing d = 100A they 
found r = 1.5 x 10~ 13 kgs _1 . An applied voltage of 1 V resulted in a friction 
r = 3 x 10~ 12 kg/s at 300 K with d = 20nm. 

In Ref. |3] the non-contact friction has also measured for fused silica sam- 
ples. Near the silica surface the friction was found to be an order of magni- 
tude larger than for the gold sample. The silica sample had been irradiated 
with 7 rays which produce E' centers (Si dangling bonds) at a density of 
7 x 10 17 cm -3 . Although the sample is electrically neutral overall, the E' cen- 
ters are known to be positively charged, creating enhanced field inhomogene- 
ity and causing the non-contact friction to rise another order of magnitude. 

Attempts to explain the observed friction in terms of the van der Waals 
friction have not met with much success since the van der Waals friction for 
good conductors like copper has been shown UHl HZj to be many orders 
of magnitude smaller than the friction observed by Stipe et.al.. In ^H] it was 
proposed that the van der Waals friction may be strongly enhanced between 
a high resistivity mica substrate and silica tip. However in |3] the mica 
substrate and silica tip were coated by gold films thick enough to completely 
screen the electrodynamic interaction between the underlying dielectrics. 

At small separation d ~ lnm, resonant photon tunneling between adsor- 
bate vibrational modes on the tip and the sample may increase the friction 
by seven order of magnitude in comparison with the good conductors sur- 
faces [injl2DI- However, the distance dependence (~ 1/d 6 ) is stronger than 
observed experimentally [3]. 

Recently, a theory of noncontact friction was suggested where the friction 
arises from Ohmic losses associated with the electromagnetic field created by 
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moving charges induced by the bias voltage |2T] • In the case of a spherical tip 
this theory predict the same weak distance dependence of the friction as ob- 
served in the experiment, but the magnitude of the friction is many orders of 
magnitude smaller than found experimentally. However, we have shown that 
the electrostatic friction can be greatly enhanced if there is an incommensu- 
rate adsorbed layer exhibiting acoustic vibrations j221- This theory gives an 
explanation for the experimentally observed bias voltage contribution to the 
non-contact friction. 

In this article we extend the theory presented in [22J to include the con- 
tribution to friction from the spatial variation of the surface potential and 
from the spatial fluctuation of the electric charge of charged defects in the 
bulk of the dielectric. We also show that the electrostatic friction as well as 
the van der Waals friction can be greatly enhanced for 2D-systems, e.g. a 
2D-electron system or an incommensurate layer of adsorbed ions exhibiting 
acoustic vibrations. The origin of this enhancement is related to the fact 
that the screening in 2D-systems is much less effective than for 3D-systems. 
An atomic force microscope tip charged by the bias voltage, or by the spa- 
tial variation of the surface potential, and moving close to the metal surface 
will induce "image" charge in the 2D-system. Because of the finite response 
time this "image" charge will lag behind the tip, and this effect result in 
force acting on the tip, referred to as the "electrostatic friction". However, 
the weaker screening effect in the 2D-system will result in a much weaker 
restoring force, which occurs when the "image charge" is displaced from the 
equilibrium position, and this result in larger lag of the "image" charge in 
2D-systems in comparison with 3D-systems. 

Another contribution to the friction from the electric field, is associated 
with the time-dependent stress acting on the surface of the surface due to 
the tip oscillations. This stress can excite acoustic phonons, or induce non- 
adiabatic time-dependent deformation. In this article we develop theories 
of phonon and internal friction due to the time- dependent stress acting on 
the surface. We show that this stress depends on the bias voltage as V 2 
resulting in to the friction coefficient r ~ V 4 . Thus this mechanism can be 
ruled out as an explanation of the experimental data observed in j^j, where 
T ~ V 2 . In the case of phonon friction only phonons with q < uj/c s can be 
excited, where q is the component of the wave-vector parallel to the surface 
of the substrate, u is the frequency of the tip oscillations, and c s is the 
sound velocity. Thus in the phase space the area occupied by the excited 
phonons ~ (uj/c s ) 2 . For electromagnetic mechanisms of the friction (which 
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include the electrostatic and van der Waals friction) all components of the 
electromagnetic field with q < d{ , where d\ is the radius of interaction (see 
below, typically d\ ~ lOOnm), are important. Thus for the metal substrate 
in the typical case (ujdi/c s ) 2 <C 1, the phonon friction is negligible in the 
comparison with the electromagnetic friction. 

2 Electrostatic friction due to a bias voltage 
and the spatial variation of the surface po- 
tential 

2.1 A general theory 

We begin by considering a model in which the tip of a metallic cantilever 
of length L is a section of a cylindrical surface with the radius of curvature 
R (Fig.l). The cantilever is perpendicular to a flat sample surface, which 
occupies the xy plane, with the z- axis pointing outside the sample. The 
tip displacement u(t) = xu e" 4w ' is assumed to be parallel to the surface 
(along the x axis), which will be a good approximation when the oscillation 
amplitudes Uq is sufficiently small. The cantilever width w, i.e. the size in the 
direction perpendicular to the xz plane, is taken to be much larger than the 
thickness c {w ^> c), and d is the separation between the tip and the sample 
surface. It is straightforward to obtain the static electric field distribution in 
the practically important case os small distances d such that the electrostatic 
field of the entire cylinder effectively the same as that due to its bottom part. 
( The criterion that d ut satisfy for this to be the case is given by \J d/ R -C 1 .) 
The problem is then reduced to solving the two-dimension Laplace equation 
with the boundary conditions that the potential has constant values V and 
at the metallic surfaces of the tip and the substrate. The electric field 
distribution outside the conductors is equal to the field due to two charged 
wires passing through points at z = ± = ±J (d + R) 2 — R? [221 • The wires 
have charges ±Q per unit length, Q = CV, where C^ 1 = 2\n[(d+R+d%)/R]. 
The electric potential at a point r exterior to the tip and sample is given by 

<^ (r) = — 2Q [In |r — r + | — In |r — r_|] 

(1) 
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where r± = ±zd\. The attractive cantilever-surface force can be calculated 
straightforwardly using Eq.(JTJ) [2T] . 

A somewhat different picture applies in the case of a oscillating charged 
tip. The cantilever charge is not changed when its tip moves parallel to the 
surface, while the sample charge varies in time at any fixed point. The electric 
field from the oscillating tip will induce an electric charge in the sample and 
this will result in to induced electric field outside the sample. The oscillating 
electric potential due to the tip oscillation at a point r exterior to the tip 
and sample is given by 

y?i(r, t) =v9i(r)e"'^ + c.c, (2) 

where 

<p x (r) =iQu Q r d rre iqx [e- |9l|2 - 2+l - e-^-^RJq.u)] , (3) 
J-oo \q\ L J 

and R p (q, u) is the reflection amplitude for the p— polarized electromagnetic 
waves. The electric field is given by E(r) = — Vyj(r). The energy dissipation 
per unit time induced by the electromagnetic field inside of the metallic 
substrate is determined by integrating the Poynting vector over the surface 
of the metal, and is given by 



£ J dSz-[E(r)xB*(r)] z=+0 + c.c. = ~ J dS ^i(r)-^(r 

AujQ 2 \u \ 2 w / dqqe' 2qdl ImRp(u,q) 
Jo 




Taking into account that the energy dissipation per unit time must be equal 
to 2uj 2 T |uo| 2 , using (@J) gives the friction coefficient: 

r = hm 2C 2 V 2 w r dgge~^ lmR ^ q \ (5) 

u;^0 Jo UJ 

Without derivation Eq.(J3]}- was firstly presented in (221 • Now we assume that 
the electric potential on the surface of the tip is inhomogeneous, consists of 
the domains or "patches". Thus the cylinder with linear size w is divided 
on smaller cylinders with the linear size Wf w = J2i w i ^ w i ^ \fdR, and 
with the surface potential Vi S = V + Vj, where V is the bias voltage and V* 
is the randomly fluctuating surface potential for the domain i. In the case 
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of a cylindrical tip geometry all domains give independent contribution to 
friction which can be obtained from Eq.(|5|) after replacement V — > V + V,i 
and w —>■ Wi . The contribution to friction from all domains is given by 

r = E r« = E ft* 2C 2 (y + vtf Wi [°° d qqe - 2 ^ lmRp{u} > q) 

= lim 2C 2 (V 2 + V 2 )w f°° dqqe- 2 ^ lmRp{uJ,<l) (6) 

where we take into account that the average value of the fluctuating surface 
potential (Vj) = J2i"WiVi = and V 2 = J2iWiV 2 /w, so that Vo is the root 
mean square variation of the surface potential. According to Eq.fjHJ), bias 
voltage and patch contributions to the friction have the same dependence 
on d. Sukenik et al, studied the root mean square variation of the surface 
potential due to thermally evaporated gold using the Stark effect in sodium 
atoms [21]. The films were partially optically transparent with a thickness 
of 42 nm and heated at 120 °C for several hours in vacuum. They deduced 
the magnitude of the fluctuating surface potential to be V = 0.15V, and 
showed that the scale of the lateral variation of the surface potential is of 
the order of the film thickness. The measurement of the non-contact friction 
between a gold tip and the gold sample gave Vo ~ 0.2V 3J thus confirming 
the prediction of the theory that this parameter is determined by the root 
mean square variation of the surface potential. 

Now, let us consider spherical tip (radius R) with the constant voltage 
surface domains with the linear size R{. If R 3> Ri ^> V dR the domain on 
the apex of the tip will give the main contribution to the friction. In this 
case we can neglect the spatial variation of the surface potential and the 
electric field induced by the bias voltage is approximately the same as that 
which would be produced in the vacuum region between two point charges 
±Q t = ±C(V + Vi) located at 



z = ±dx = ±\/3Rd/2 + sJ{?>Rd/2) 2 + Rd 3 + d A (7) 
where 

C^-^ (8) 

2d K ' 

It can been shown that the electrostatic force between the tip and the metal 
surface within this approximation agrees very well with the exact expres- 
sion for a sphere above a metal surface (2H1- The vibrations of the tip will 
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produce an oscillating electromagnetic field, which in the vacuum region co- 
incides with the electromagnetic field of an oscillating point charge. The 
friction coefficient for a point charge moving parallel to the surface due to 
the electromagnetic energy losses inside the sample, is determined by [2E] 

r^iim^f^V^^MM) (9) 

11 «->o 2 Jo to V ' 

For motion normal to the surface, T± = 2Tn. Thus, just as for the cylindrical 
tip geometry, for a spherical tip the friction depends parabolically on the 
bias voltage. However for a spherical tip the parabola begins from zero in 
contrast to a cylindrical tip, where the parabola begins from a finite positive 
value. 



2.2 Clean surface 

For a clean flat surfaces the reflection coefficient is determined by the well- 
known Fresnel formula 

R P = (io) 

In this case, for the tip radius R 3> d and for a metal with the dielectric 
function e = 1 + inia/uj, where a is the conductivity, Eq.© gives: 

Neglecting the contribution from the spatial variation of the surface potential, 
this formula was obtained recently in [21] using a less general approach. With 
w = 7 ■ 10~ 6 m and a = 4 • 10 17 s _1 ( corresponds to gold at 300 K), and with 
d = 20nm and V = lVolt, Eq.dTT)) gives T = 2.4 • 10~ 20 kg/s which is eight 
orders of magnitude smaller than the experimental value 3 • 10~ 12 kg/s |3]. 

Assuming R » causing Q and (fTUj) gives the friction between a spherical 
tip and a clean sample surface 

ol/2 D1/2T/2 

_ 3 R V 

This expression is only a factor 1.6 smaller that the result obtained inde- 
pendently in [21] using a less general approach. For the same parameters 
as above and at d — 20nm, the friction for a spherical tip is two order of 
magnitude smaller than for the cylindrical tip. 
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To get insight into possible mechanisms of the enhancement of non- 
contact friction it is instructive to note that qualitatively Eq.([ll|) can be 
obtained from the following simple geometrical arguments [27J. The vibrat- 
ing tip will induce current in the sample in a volume with the spatial dimen- 
sions L x , L y and L z . The instantaneous dissipated power in the sample is 
given by P ~ I 2 r, where / is the current and r is the effective resistivity. 
The current / is proportional to the tip velocity v x , and can be written as 
I ~ v x Q t /L x , where Q t is the charge of the tip. The effective resistance r 
can be approximated by the macroscopic relation r = pL x /L y L z , where p 
is the resistivity. Using this simple expressions for current / and resistance, 
and using the relation Q t = C t V s (where C t is the tip-sample capacitance) 
for the induced charge, the instantaneous power dissipation is 

? ,2^2t/2 

P = I 2 r ~ p-fff- (13) 

Comparing this expression with P = Tv 2 we get 

C 2 V 2 

r ~ pjYt ^ 

Ll X Llyij Z 

For a cylindrical tip vibrating above the clean surface L y ~ w and L x ~ L z ~ 



d\. For d R the tip-sample capacitance C t ~ wyR/8d and d\ ~ \ / 2dR. 
Substituting these expressions in Eq.(^3J) gives Eq.fjlljl to within a numerical 
factor of order of unity. From Eq. (fUjl it follows that the friction will increase 
when the thickness L z of "dissipation volume" decreases. This is the reason 
for why 2D-systems may exhibit higher friction than 3D-systems. 



2.3 Film on-top of a high-resistivity substrate 

^From the qualitative arguments given above it follows that for a thin metal 
film on-top of a high resistivity substrate, e.g. a dielectric or a high resistivity 
metal, the friction will be larger, than for an infinitely thick film. In this 
case the thickness L z of the volume, where the dissipation occurs, will be 
determined by the thickness of the film, and according to Eq. (jl4j) this will 
give rise to a strong enhancement of the friction. 

For a planar film with thickness df and dielectric constant e 2 on-top of a 
substrate with dielectric constant e 3 , the reflection coefficient is determined 



9 



by 

R = R P 2i - R P 23exp(-2qd f ) 
p 1 - Rp2iRp23 exp(-2qd f ) 

where 

R pij = 7~7 i ( 16 ) 
tj -t- tj 

where index 1 is associated with vacuum. For a metallic film on a dielectric 
substrate, or a metallic film on a metallic substrate with o<i ^> a 3 , for g?i 3> (i/ 
and R^$> d Eqs. © and (fT5|) gives 

w (F 2 + y 2 )it!V 2 

/ 2 9 /%ac//c/ 3 / 2 ' 1 J 

This is greater by a factor of 2d\/df than the corresponding friction for the 
infinitely thick sample. For thin film the effective resistivity of the substrate 
is increased, giving rise to additional ohmic dissipation. In [21] Eq. (jl7j) was 
obtained using a less general approach and neglecting the spatial variation 
of the surface potential. The conditions necessary for the validity of Eq. (jl7j) 
could not be determined in this simplified approach. 

2.4 2D-system on-top of a dielectric or metal substrate. 

Let us now consider a 2D-system, e.g. electronic surface states or a quantum 
well, or an incommensurate layer of ions adsorbed on a metal surface. For 
example, for the Cs/Cu(100) system experiment suggests the existence of an 
acoustic film mode even for the very dilute phase (6 ~ 0.1). This implies 
that the Cs/Cu(100) adsorbate layer experience a negligible surface pinning 
potential. The reflection coefficient for p-polarized electromagnetic waves 
can be obtained using the approach proposed in |2E1- This gives (see 
detailed derivation is given in Appendix A) : 

_ 1 - 1/e + Airqn a a\\/e - qa(l - 47in a qa\\) 
1 + 1/e + Anqn a aii/e + qa(l + ATm a qa\\) ' 

where n a is the concentration of the free carries of the charge per unit area. 
The polarizability an for the 2D-system in the direction parallel to the surface 
is taken to be 

e* 2 

a " = ~M(u; 2 + ^||)' (19) 
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where 7711 is the damping constant, e* and M are the effective charge and the 
mass of the moving particles, respectively. In comparison with the expression 
obtained in [2H], Eq. (jl8|) takes into account that the 2D-system is located a 
distance a away from the image plane of the metal. Although this correction 
to the reflection coefficient is of order qa <C 1 , for a 2D-system on-top of a 
good conductors ( |e| ^> 1), it gives the most important contribution to the 
energy dissipation. 

For good metals (|e| ^> 1), from Eq.ffTHj) we get 

Im.R « — — ^ q , (20) 

[U} z — UJ Z J Z + UJ Z 7] Z 

where u 2 = 47m a e* 2 aq 2 /M. In the case of a 2D-structure on-top of a dielec- 
tric, the factor qa in Eq.([2U|) and in the expression for uj 2 must be replaced 
by 1/e, where e is the dielectric function of the substrate. Using (j2*U|l in © 
for R 3> d we get 

_ wr ] MRi/ 2 (V 2 + V 2 ) 
ad 2 9 / 2 d 3 / 2 7m a e* 2 ' 1 ' 

This friction exhibits the same distance dependence as observed experimen- 
tally [3 . The same expression for the friction is valid for a 2D- structure 
on-top of a dielectric. Comparing Eqs. ljllj) and (|2"T|) we find that a 2D- struc- 
ture on-top of a substrate gives the same friction as for the clean surface with 
the effective conductivity a e ff = n a e* 2 jMr(ld\. We obtain agreement with 
experiment at d — 20nm if a e ff m 4 • 10 9 s -1 . In the case of a 2D-electron 
system, for R = 1/im such an effective conductivity is obtained if rj = 10 14 s _1 
and n a = 10 15 m~ 2 . For Cs/Cu(100), for n a = 10 18 m" 2 (9 w 0.1) the electric 
charge of the Cs ions e* = 0.28e Due to the similarities of Cu and Au 
surfaces, a similar effective charge can be expected for the Cs/Au surface. For 
such a 2D- system agreement with experiment is obtained for n a = 10 18 m -2 
and 77 = 10 n s _1 . In [22] we estimated the damping parameter for a Cs atom 
associated with the covalent bond r]\\ cov = 3 ■ 10 9 s _1 However the colli- 
sions between the ions, and between the ions and other surface defects, will 



also contribute to rj. In this case r] co i ~ vt/1 where vt ~ y^sT/M, and I is 
the ion mean free path. For T = 293K and / ~ lnm we get r) co i = 10 11 s _1 . 

For a spherical tip, with a 2D-system on-top of the substrate, from 
Eqs.(J2nj) and © for R ^> d we get the contribution to the friction from 
the 2D-system 

3RM.V 2 

ad 2 6 d7rn a e* 2 ' y } 



11 



At d = 20nm this friction is ~ two order of magnitude smaller than for the 
cylindrical tip. 

2.5 Friction due to spatial fluctuations of static charge 
in the bulk of the sample 

In this section we consider a dielectric substrate with a stationary, inhomo- 
geneous distribution of charged defects. Such a situation was investigated 
experimentally |3| by employing a fused silica sample irradiated with 7 rays. 
In the course of irradiation, positively charged centers (Si dangling bonds) 
are generated. Randomly distributed positive charges are compensated by 
randomly distributed negative charges, thus on average the sample is elec- 
trically neutral. We model the sample as consisting of microscopically small 
volume elements AV^. Each element is chosen sufficiently small that not 
more than one charge center is present in it. The electric charge of each 
element is equal to ±e or 0, in such away that the average (qi) = 0. We 
will consider the fluctuations of charges in different volume element i,j to 
be statistically independent, so that ((Mj) = for i 7^ j. The mean square 
of charge fluctuations within a given element ((Mi) ~ 2ne 2 , where n is the 
average number of positive charges in one volume element. In the absence 
of the cross terms the average tip-sample friction coefficient is determined 
by adding friction coefficient from all charges qi. According to Eq.Q, the 
contribution to the friction coefficient from charge qi in the element AVi is 
given by 

ATjii = lim ne 2 [°° dqq 2 e~ 2 ^ lmR P^^^ (23 ) 
^->o Jo uj 

where di = D(xi,yi) — z%. Here the coordinates Xi,yi,Zi give the position 
of the i-th volume element in the substrate, and D(xi,yi) is the distance 
between the substrate and points Xi,yi located on the surfaces of the tip. 
The total friction coefficient is obtained by summing over all the elements. 
Replacing the sum by an integral (nJ2 — * cj d 3 r, where c is the number of 
the positive charge centers per unit volume), and integration over z gives 

r„ = hm ^ r dqq [ dx [ (24) 
<^o 2 Jo J J LU 

For a cylindrical tip D(x, y) = d + x 2 /2R, and we get 
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q = lim jT d qq y^ Im ^ q) (25) 

Using the same parameters as in Sec. 2. 2, for a gold tip separated by d = 
lOnm from a dielectric sample with c = 7 ■ 10 17 cm -3 we get Tii = 4.4- 10~ 20 kg 



s- 1 . 



For the tip surface with a 2D-structure on it, using Eq.(j20J) we get 

„ 1 / e \ 2 R cw „ , e 2 cw , , 

r -« = ^( ? ) V^ M " = T6^ (26) 

With a eff = n a e* 2 /2Mr]d 1 = 4 • 10 9 , c = 7 • 10 17 cm~ 3 , and with the other 
parameters the same as before, we get for d = lOnm, r^y^ 3.5-10~ 12 kg s -1 , 
which is nearly the same as was observed experimentally pH . Thus our theory 
of friction between a gold tip and silica substrate with an inhomogeneous 
distribution of the charged defects is consistent with the theory of friction 
between a gold tip and gold substrate ( see Section 2.4). In both theories we 
have assumed that the gold surfaces are covered by a 2D-structure. 

The study above has ignored the screening of the electric field in the 
dielectric substrate. This can be justified in the case of very small tip- sample 
separations (substantially smaller than screening length), as only defects in 
the surface layer of thickness d contribute to the integral in Eq. ()24|) . When 
the screening is important, the effective electric field outside the sample will 
be decreased by the factor (e + l)/2 [23], and the friction coefficient will be 
decreased by the factor ((e + l)/2) 2 , which is equal to ^6.25 in the case of 
silica. However, the inhomogeneity of the surface of the tip may be larger 
than that of the sample surface, so that the damping parameter 77 may be 
larger for the 2D-structure on the surface of the tip. This increase in rj and 
screening effects will compensate each other. 



3 Van der Waals friction 

In this section we consider the van der Waals friction between two surfaces 
covered by 2D-systems. The frictional stress between two flat surfaces to 
linear order in the relative velocity v can be written in the form: a = jv. 
According to in the case of the van der Waals friction the contribution 
to the friction coefficient 711 from the p-polarized electromagnetic waves is 
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given by 

h r°° , / dn 



711 



2tt 2 



Jo \ ouj / jo 



x Imfi "' m ^ |i-e-^ lP fi 2 pi 2 (27) 

where and i?2 P are the reflections coefficients for the surfaces, and n = 
[exp(huj/k B T) — In [1211201 we have shown that resonant photon tun- 

neling between two Cu(100) surfaces separated by d = lnm and covered by 
a low concentration of potassium atoms gives rise to a friction six orders 
of the magnitude larger than for clean surfaces. The adsorbate induced en- 
hancement of the van der Waals friction is even larger for Cs adsorption on 
Cu(100). In this case, even at low Cs coverage {6 ~ 0.1), the adsorbed layer 
exhibit an acoustic branch for vibrations parallel to the surface [2H|, and 
according to Eq. (fTS|) . at small frequencies the reflection coefficient is given 
by 

2qau 2 

L)' Z — 0J 2 + iuJT\ 

where uj\ = 4nn a e* 2 aq 2 /M. Using Eq. (j23) in Eq. (J2ZJ) for 



flp = i- .. 2 hrr— (28) 



a 4irn a e* 2 a 

Md 2 ^ ' 



gives 



kuTa 2 , , 

7|l « 0.62^. (29) 

It is interesting to note that according to 7y does not depend on n a , e*, 
and M. However, Eq. (J28|) is only valid when there are acoustic vibrations 
in the adsorbed layer. For Cs adsorbed on Cu(100) the acoustic vibrations 
exist only for 9 > 0.1 29J. The friction coefficient for a cylindrical atomic 
force microscope tip can be estimated using (HOI EE] 

f°° , / / \\ k B Ta 2 R°- 5 w 
2w ' ' " 1 "'•><- 



./ o "fc 7 ,(,M) = 0.68 ^°f 5 W (30) 



where i? is the radius of the curvature of the tip and w is its width, and 
^y\\(z(x)) the friction coefficient between two flat surfaces at the separation 
z(x) — d + x 2 /2R. In Section 2 we have shown that the experimental data 
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in [3] can be explained by assuming that the gold surfaces are covered by 
adsorbed layer of ions like Cs on Cu(100) with the damping constant rj « 
10 11 s~ 1 . With this value of rj and using a = 2.94A R — 1/im, w = 
7/xm, T = 293 K we find that if d < 3nm the contribution from the van 
der Waals friction will dominate over the contribution from the electrostatic 
friction. However, in the experiment a strong enhancement in the friction 
was not observed at such short separation. Thus, most likely a 2D-system 
of electronic origin is responsible for the enhancement of the electrostatic 
friction. In this case (see Section 2) rj e i ~ 10 14 s _1 and the van der Waals 
friction will give a negligible contribution for practically all separations. Fig. 2 
shows how the friction between the copper tip and the copper substrate 
depends on the distance d , when the surfaces of the tip and the substrate 
are covered by a low concentration of the Cs atoms, and for clean surfaces. 
In comparison, the friction between two clean surfaces at the separation d = 
lnm is eleven orders of the magnitude smaller. However, the friction between 
clean surfaces shown on Fig. 2 was calculated in the local optic approximation. 
For parallel relative motion non-local optic effects are very important [20 , 
and when it is taken into account, at d — lnm the friction between adsorbate 
covered surfaces will be seven orders of the magnitude larger than the friction 
between clean surfaces. 



4 Phonon and internal non-contact friction 

4.1 Non-contact friction due to excitation of substrate 
phonons 

Consider a tip which performs harmonic oscillation, u = uq exp(— iut) + c.c, 
above an elastic body with a flat surface. This will results in a fluctuating 
stress acting on the surface of the solid which excite acoustic waves with 
parallel wave number q < u>/c s , where c s is the sound velocity. The stress 
o iz acting on the surface of the elastic solid can be represented through the 
Fourier integral 

72 

Mx, t) = j —^a^uoe^-^ + c.c. (31) 

Using the theory of elasticity (assuming an isotropic elastic medium for sim- 
plicity), one can calculate the displacement field Ui on the surface z = in 
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response to the surface stress distribution a iz 

«i(x, t)= j jA^M^q, ufa (q)n e iqx -^ + c.c. (32) 

The energy dissipation per unit time equals 

P = J d 2 x(Ui (x,t)cr i2 ((x,t)) = 

2iU I J^ l ™ M *M^)<^fa^)Wo\ 2 (33) 

where (...) stands for the time averaging. The explicit form of the stress 
tensor in the model of the elastic continuum is given in [22] (see also Appendix 
B). The energy dissipation per unit time must be equal to V <ii (t) 2 >= 
T2u 2 \uq\ 2 . Comparing of this expression with ()33|) gives 



r d 2 q ImMy(q,w) 
T = Jj2^ u ^ (q) ^ (q) (34) 

At typical experimental conditions we have u ~ 10 3 — 10 6 s _1 and qr* < 
ur* fc s < 1CT 3 << 1, where effective radius of the interaction r* m y/dR, and 
where d is the separation between the tip and the sample, and R is the radius 
of curvature of the tip. In Appendix B it was shown that in this case the con- 
tribution to the friction from excitation of acoustic waves can be determined 
by calculating the energy dissipation due to oscillating point force applied to 
the surface of the semi-infinite elastic continuum. These calculations were 
done in the connection with the vibrational energy relaxation of adsorbates 
|33j . According to this theory the friction coefficient for vibration of the tip 
normal to the surface is given by 

UK 2 
An pel 

where m 1.65, c t is the transverse sound velocity of the solid, p is the mass 
density of the sample, K = dF/dd, where F(d) is the force acting on the tip 
due to interaction with the sample. 

In Appendix B it was shown that for vibration of the tip parallel to the 
flat surface the friction coefficient due to excitation of the acoustic waves is 
given by 

47T pCl 



^F 2 M (36) 
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where £n ~ 1.50. From the comparison of the Eqs.(|35l- we get that 
rn/r_L ~ (uod/c t ) 2 <C 1. We consider now two different contributions to the 
tip-sample interaction. 

4.1.1 Van der Waals interaction 

Accordingly to the Lifshitz theory [Hj the stress a zz (d) acting on the surface 
of two identical semi- infinite bodies due to van der Waals interaction at 
small separation d <C c/u p (where u p is the plasma frequency) and d <C At 
is given by: 

^ = I^f d #F$- (37) 

In the Drude model the explicit form of e is 

«(« = 1 + (38) 

For typical metal the damping constant rj <C uj p and can be neglected when 
integrating Eq.(pT7j). It follows from Eqs. (j3*7j) and (jSHj) that 

a zz = -Jp— (39) 

For the spherical tip of radius R using the same approximation as in Eq. (j30|) 

we get 

™ = S < 4 °> 

and 

K' = ™p- (41) 
Similarly, in the case of a cylindrical tip we have 

and 

K = V<F,2 ( 43 ) 
For copper tip separated from a copper substrate by d = lOnm, and with 
R = 1/xm, w = 7/xm, we get for spherical tip = 6.3 ■ 10 _18 kg s _1 and for 
cylindrical tip = 1.3 • 10~ 14 kgs -1 . The phononic friction decreases as d~ e 
and d~ 7 for spherical and cylindrical tip, respectively. 
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4.1.2 Electrostatic interaction due to a bias voltage 



In the presence of the bias voltage V the attractive force between the tip and 
the sample at d <C R is given by 

nd) = (44) 

for a cylindrical tip, and 

DT/2 

n*> - (45) 

for a spherical tip. For bias voltage V = lVolt, and with the other parameters 
the same as above, we get T s ± = 8.8 • 10~ 17 kg s _1 and T c ± = 1.2 • 10~ 13 kgs _1 
for the spherical and cylindrical tip, respectively. Note that in this case the 
friction depends on the bias voltage as V A . 

For the vibrations of the tip parallel to the sample surface the expression 
for the friction coefficient contains the addition small factor (u)d/c s ) 2 <C 1. 
Thus the friction coefficient for parallel vibrations of the tip will be by many 
orders of magnitude smaller than for normal vibrations. 



4.2 Non-contact friction due to internal friction of the 
substrate 

In studying of the phononic friction in Section 4.1 it was assumed that the 
deformations of the solids are purely elastic. However, the deformation will 
be purely elastic or adiabatic only for infinitesimally small velocity, so that 
at every moment of the time the system stays in the equilibrium state. How- 
ever, real motion always occurs with finite velocity, and the body does not 
stay in equilibrium; and thus "flow-processes" occur, which tend to bring it 
back to equilibrium. This leads to non-adiabatic deformations, resulting in 
dissipation of the mechanical energy. 

The energy dissipation is determined by two kind of processes. First, 
in the presence of a temperature gradient in the body, result in heat flow. 
Secondly, if in the body occurs some kind of internal motion, than non- 
adiabatic processes occur, related with finite velocity of the motion; these 
processes of energy dissipation can be denoted, as in liquids, as internal 
friction or viscosity. 

The friction coefficient due to the internal friction is determined by Eq. (j34j) . 
However, in contrast to the phononic friction, large values of q ^> u/ct play 
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the most important role for the internal friction. For q ^> uo/c s the tensor 
component M zz is given by [32J 

M„ = ^ (46) 

where E{uS) is the complex elastic modulus and v is the Poisson ratio. 

4.2.1 Van der Waals interaction 

For R 3> d only the a zz component of the stress tensor due to the van der 
Waals interaction is important. In this case, for vibrations of the cylindrical 
tip parallel to the sample surface, we get 

a z (q) = J d 2 xe l ^—a 



zz \ 



= _ ^ B in(^/2) (3 + e2 + 30e _ s (47) 

where £ = \/2dRq x . Using (j4"Tj) and (j4T?j) in (J34"j) we get for a cylindrical tip 

_ 75tt wh 2 u 2 p Im(£/(1 - v 2 )) 



11 2 16 d e u\E/(l-u 2 )\< 
For the spherical tip similar calculations give 



(48) 



_ 0.25 R^hW^miE/jl-u 2 )) 
11 2 9 V2n d 11 ' 2 oo\E/(l~u 2 )\ 2 ' 1 ' 

In general, Im[£'(u;)/(1 — v 2 )\ has many resonance peaks, corresponding to 
different thermally activated relaxation processes. One important source of 
internal friction at high frequencies is related to thermal currents: elastic 
compression of a material is commonly associated with heating effects. If the 
compression takes place sufficiently rapidly, there is no opportunity for heat 
to be conducted away, while for very slow compression temperature gradi- 
ents are eliminated by thermal conduction. In both these cases the process of 
compression will be reversible. In the former case it will be adiabatic and in 
the latter one - isothermal. In both these limiting cases the contribution from 
thermal current to the internal friction will be negligible. However, in the 
intermediate frequency regime we expect dissipation of mechanical energy 
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into heat. The characteristic frequency for the maximum dissipation will be 
of order uo t = 1/r where, from dimensional arguments, we expect the relax- 
ation time r ~ l 2 /D, where I is the linear size of the compression region and 
D the thermal diffusibility D = k/ ' pC v (where C p is the specific heat and k 
the heat conductivity). For I ~ 10 3 A, this gives for gold Ut ~ 10 u s~ 1 , which 
is much higher than the resonance frequency of the cantilever of the atomic 
force microscope. Another very important contribution to the internal fric- 
tion is point-defect flipping. This involves thermally activated transitions of 
point defects or loose sites in crystalline and amorphous network. A special 
case is the vibrational motion of adsorbates on the surface of the substrate 
and/or on the tip, as was treated separately above. Another contribution 
to the internal friction comes from grain-boundary slip [31] . For a copper 
cylindrical tip and a copper substrate using d = lOnm, w = 7/zm, R = 1/im, 
uj = 10 4 s _1 , and, as is typical for metals [35], lm.E(u)/ \E(u)\ ~ 10~ 5 and 
E pa 10 11 N/m 2 , gives Tfj ph 10 _16 kg-s _1 . Thus at this separation the inter- 
nal friction gives much smaller contribution to the friction coefficient than 
electrostatic friction due to bias voltage or spatial variation of the surface 
potential. However, internal friction can give the dominant contribution for 
small separation d < lnm. For the spherical tip with R = 1/im the friction 
coefficient is two order of the magnitude smaller. Finally we note, as a cu- 
riosity, that the internal friction of solids gives a very important contribution 
to the rolling resistance of the most solids [3B], and is the main contribu- 
tion to rubber friction on rough substrates, e.g. road surface (HE], where in 
the transition region between the rubbery and glassy region of the rubber 
visco-elastic spectra, \mE{u)/ \E(u)\ ~ 1. 



5 Summary 

We have studied how the electrostatic friction between an atomic force mi- 
croscope tip and a substrate depends on: (a) the bias voltage, (b) the spatial 
variation of the surface potential, and (c) the spatial fluctuation of electric 
charge. We have found that the electrostatic friction can be greatly enhanced 
in presence of a 2D-system on the surface of the sample or on the tip. On 
metal surfaces such 2D-system can result from surface electronic states, or 
from an incommensurate layer of adsorbed ions. We have shown that the ex- 
perimental data observed in [3] can be explained by the electrostatic friction 
in presence of such a 2D-system. The theory predicts the same magnitude, 
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distance and bias voltage dependence of the friction coefficient as it was 
observed in the experiment [3], and explains the bias- voltage-independent 
contribution to friction. The theory of friction between a gold tip and sil- 
ica substrate with an inhomogeneous distribution of the charged defects is 
consistent with the theory of friction between a gold tip and gold substrate. 
In both theories we have assumed that the gold surfaces are covered by 2D- 
structure. 

The electrostatic friction was compared with the van der Waals friction 
arising from quantum and thermal fluctuations of the current densities inside 
the bodies. The van der Waals friction as well as the electrostatic friction 
can be greatly enhanced in presence of identical 2D-system on the surfaces 
of the tip and the substrate. The van der Waals friction is characterized by a 
stronger distance dependence than the electrostatic friction, and may domi- 
nate at small separation. The van der Waals friction between 2D-systems can 
be so large that it can be measured with present state-of-the-art equipment. 

Phonon and internal friction can be ruled out as mechanisms responsible 
for non-contact friction observed in [Sj because they predict stronger distance 
and bias voltage dependence. For metal substrate the phonon friction asso- 
ciated with excitation of acoustic phonons is negligibly small in comparison 
with the electromagnetic friction (especially for motion of the tip parallel to 
the substrate surface) because of small area in the phase space available for 
these phonons. 
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A Fresnel reflectivity for ^-polarized electro- 
magnetic waves with 2D-structure correc- 
tions 

We consider a semi-infinite metal having a flat surface which coincides with 
the xy plane, and with the z axis pointed along the inward normal. The 
metal surface is covered by an adsorbate layer located at z — —a. Let the 
xz plane be the plane of incidence of evanescent electromagnetic plane wave, 
with the parallel component of the wave vector q pointed along the x-axis. 
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The macroscopic electric field takes the form 

Ie~ pz + Ke pz 
E = ( x { Ae- pz + Be pz 

Te -sz 



z < —a 
-a < z < 
z > 



(50) 



where p = (q 2 — (cu/c) 2 ) 1 ^ 2 , s = (q 2 — (u / ' c) 2 s{uj)) 1 ^ 2 , and e is the dielectric 
function of the metal. According to [2E] the boundary conditions at z = —a 
can be written in the form 



B 7 e 



A z e pa + B~ pa - R z e~ pa - I z e pa = Anpn a a\\ (R z e~ pa - I z e pa ) (51) 

pa_ Ae pa_ R - P a + I P a _ _47W*± 



P 



[R z e~ pa + I z e pa ) (52) 



where anm is the polarizability of the adsorbate in the direction parallel 
(normal) to the surface. From the ordinary boundary conditions at z — it 
follow 

B z = -!—-^-A z (53) 
ep + s 

For a 2D-system a±_ = 0, and for q ^> uj/c and qa <C 1 Eqs. (j51H5^|) give the 
reflection coefficient R z , given by Eq.( ITHj). 



B Friction coefficient due to excitation of the 
acoustic waves 

According to |S2] the tensor M in Eq.(|H^|) is given by 

i / 1 



pet K S(q,u 



10 



-h| — ] (pizz + p t 



Q(q,u)(zq_-qz 
1 



+ nn- 



Pt 



where q = q/g, n = z x q, and where 



5 

Pt 



\^j-2q 2 j +4q 2 p t pi, 
2q 2 - uj 2 /c 2 + 2p tPl , 



q 2 , pi 



(54) 

(55) 
(56) 

(57) 
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In the equations above, p, c t} and q are the mass density and the transverse 
and longitudinal sound velocities of the solid, respectively. Note that q and 
c\ are in general complex frequency dependent quantities given by 

C < = 2^)' (58) 
c 2 - E{1 ~ U) (59) 

where E(u) is the complex elastic modulus and v is the Poisson ration. 

The acoustic waves have wave number q < uj/c t . At typical experimental 
condition the frequency of the vibrations of the tip u ~ 10 3 — 10 6 s _1 and 
qRint < ^Rint/ct <C 1, where Ri nt ~ \fd~R is the radius of the interaction 
of the tip with the sample surface. In this case for the vibrations of the tip 
normal to the surface we get 

5lZ I d2X dd a ° zzi ^ d) = dd FM (60) 
where of 2 is the static stress acting on the surface of the sample. Using 
Eqs.flHJ) and ^ in Eq.flMJ) we get 



An pcf 



(61) 



where £j_ = £_u + £_u + £± s , K = dF z /dd and where the contributions from 
the longitudinal the transverse £± t , and surface (Rayleigh) acoustic 
waves are given by 



r t/ci , yW-f fR0 , 

€±l = / dx : — (62) 

(1 - 2xf + Ax J 1 - x)J(ct Cl ) 2 - x 







f 1 Ax[x-(c t / Cl ) 2 }VT^ 
Ut ^ h to (1 - 2xY + lto*[x - (c t / Cl y](l - x) ' 

£± s = nJx~^(c~Jc~iF/f(x c ), (64) 



where 



/(ar) = 4xv / ^T V /;r - ( c */ c ') 2 - ( 2x - X ) 2 > ( 65 ) 
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and where x c is the solution of the equation f(x) = 0, /' = df(x)/dx. In 
Eqs.dnHHZni) the sound velocities c t and q are assumed real, taken at uj = 0. 

For the vibrations of the tip parallel to the surface the main contribu- 
tion to the energy dissipation due to excitation of the acoustic waves gives 
component of <7, which acts in the ^-direction. For this component we get 

a h {q) = J d 2 xe^—a° zz (*) ~ k*F z {<t) (66) 
Using Eqs.(jnSl) and in Eq-dSH) we get 

r„ = |£fJH) (67) 

where ^ = % +£\\ t + Z\\ s , 

- r-te . ( 68) 

(i _ 2x ) 2 + 4x^(1 - x) ^(c t /ci) 2 - x 

^ 4/c (1-2x) 4 + 16x 2 [x-(q/q) 2 ](1-x)' 1 ; 

£_l s = ixx c \jx c - (c t /ci) 2 /f(x c ), (70) 
For most metals c t /q ~ 1/2 and for this case £j_ = 1.62 and £n = 1.50. 

FIGURE CAPTIONS 

Fig. 1. Scheme of the tip-sample system. The tip shape is characterized 
by its length L and the tip radius of curvature R. 

Fig. 2. The friction coefficient associated with the van der Waals friction 
between a copper tip and a copper substrate, both covered by low concentra- 
tion of cesium atoms, as a function of the separation d. The cylindrical tip is 
characterized by the radius of curvature R = l/j,m and the width w = 7/im. 
The other parameters correspond to Cs adsorbed on Cu(100) at the concen- 
tration n a = 10 18 m" 2 (coverage 9 « 0.1) [2211211: e* = 0.28e, rj = 10 u s _1 , 
a = 2.94A, T = 293K. (The base of the logarithm is 10) 
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